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Abstract 

In this article we explore the R-symmetry of the (2,2) 2d Wess-Zumino 
model. We study whether or not this symmetry is approximately realized in 
the Q-exact lattice version of this theory. Our study is nonperturbative: it 
relies on Monte Carlo simulations with dynamical fermions. Irrelevant opera- 
tors in the lattice action explicitly break the R-symmetry. In spite of this, it 
is found to be a symmetry of the effective potential. We find nonperturbative 
evidence that the nonrenormalization theorem of the continuum theory is re- 
covered in the continuum limit; e.g., there is no additive mass renormalization. 
In our simulations we find that Fourier acceleration of the hybrid Monte Carlo 
algorithm allows us to avoid difficulties with critical slowing-down. 
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1 Motivation 



The continuum (2,2) 2d Wess-Zumino (2dWZ) model (obtained from a dimensional re- 
duction of the 4d Wess-Zumino model [T]) is supposed to provide a Landau- Ginzburg 
description of the minimal discrete series of A/" = 2 superconformal field theories j2] • In 
the present article, we examine an important aspect of the simplest of these models — 
the one with a cubic superpotential — in the context of class of lattice actions that 
have an exact lattice supersymmetry. These lattice actions were first formulated in 
OE] using Nicolai map {S] methods, relying on earlier Hamiltonian jH] and continuum 
[7j studies that also utilized the Nicolai map. Detailed studies of the spacetime lattice 
system were performed in |SI by stochastic quantization methods and in jS] by the 
Monte Carlo simulation approach. 

Once auxiliary fields are introduced, the lattice action takes a Q-exact form: 
S = QX, as was emphasized in topological interpretation of jTU] and the lattice 
superfield approach of jTTj. Here Q is a lattice supercharge with derivatives realized 
through discrete difference operators; with respect to a discrete approximation of 
the continuum theory superalgebra, = is a nilpotent subalgebra. Because 5* is 
Q-exact, the action is trivially invariant with respect to this lattice supersymmetry: 
QS = Q^X = 0. 

In the massive continuum theory, as will be shown below, there is an exact Z2{R) 
symmetry. It is an R- symmetry, meaning it does not commute with the supercharges. 
This symmetry is spontaneously broken at infinite volume. In the massless case, i.e., 
in the critical domain, the classical R-symmetry is enlarged to U{1)r. It cannot be 
spontaneously broken since it is a continuous symmetry in 2d [I^ • If the lattice theory 
has the correct continuum limit, it should reproduce these features. On the other 
hand, these R-symmetries are only approximate in the Q-exact lattice action; the 
symmetry is explicitly broken by the Wilson mass term that is used to lift doublers.^ 

It has been shown in that the continuum limit of the lattice perturbation 
series is identical to that of the continuum theory, due to cancellations that follow 
from = 0. Thus, the Q-exact spacetime lattice has behavior that is similar to what 
was found on the Q, Q^-preserving spatial lattice in 0. However, it was also shown 
in ^T] that the most general continuum effective action that is consistent with the 
symmetries of the bare lattice action is not the (2,2) 2d Wess-Zumino model. This 
raises the question of whether or not the good behavior of perturbation theory persists 
at a nonperturbative level. The Monte Carlo simulation results of Catterall et al. give 
hope that the desired continuum limit is obtained beyond perturbation theory [Hj. If 
so, this would be one of the few examples of a supersymmetric field theory that can 
be latticized and studied nonperturbatively by Monte Carlo simulation without the 
need for fine-tuning of counterterms. 

Though not directly related, we pause to mention that lattice super- Yang-Mills 

^This is directly related to the breaking of the so-called U{l)v symmetry, that was pointed out 

inim. 
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theories exist that do not require fine-tuning of counterterms ^H), at least within 
perturbation theory. However, in some of those cases, it is known that the fermion 
determinant has a complex phase that depends on the boson configuration T^. The 
complex phase poses a severe difficulty for Monte Carlo simulations of those sys- 
tems [15] . 

If we can show that features of the continuum theory associated with the R- 
symmetry are recovered in the continuum limit, it provides further evidence that 
the correct theory is obtained. The symmetry that we study persists in the infrared 
effective theory in a strongly coupled regime. Thus, we are testing aspects of the 
lattice theory that lie beyond the reach of perturbation theory. 

We now summarize the rest of this article: 

• In Section 121 we briefiy describe the continuum theory that the lattice model is 
supposed to define. We discuss both the general (2,2) 2d Wess-Zumino model, 
and the specific case that we select for further study: the model with a cubic 
superpotential. In that case, we show that there is a classical Z2{R) symmetry. 
We explain that it must be spontaneously broken, according to a well-known 
theorem for 2d field theories with a stable potential. 

• In Section El the lattice action is described. Symmetry aspects of the action 
are emphasized. In particular, we show that when doublers are lifted by a 
super symmetric version of the Wilson mass term, the Z2{R) is explicitly bro- 
ken. However, it is an approximate classical symmetry for the long wavelength 
modes. 

• In Section HI observables are defined that are used in our study of the critical 
domain. Basic data regarding their behavior is briefiy described. Degenerate 
minima of the effective potential lead to tunneling in the simulations at finite 
volume. We describe how observables are chosen that are insensitive to this 
finite volume effect. Next the Z2{R) symmetry is examined in our simulations. 
We find that the effective potential exhibits this symmetry to a very good ap- 
proximation. We show that as the volume is increased, spontaneous symmetry- 
breaking of Z2{R) is clearly evident. Finally, we demonstrate that the f/(l)K 
symmetry of the effective potential is recovered in the critical domain. We find 
evidence that the nonrenormalization theorems of the continuum theory hold 
to a good approximation in the limit of small lattice spacing. 

• In Section we conclude with an interpretation of our simulation results. We 
also outline future directions for research in this lattice system. 

• In Appendix the simulation methods are outlined. We describe the extent to 
which we were able to overcome critical slowing-down using Fourier acceleration 
methods, and some specifics about the algorithm in this regard. 
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2 The target theory 



We begin our discussion by considering the continuum theory. Our focus will be on 
the simplest (2,2) 2d WZ model: the one with a cubic superpotential interaction. 
It is a supersymmetric cousin to an ordinary 2d complex scalar field model. Two 
2d Majorana fermions {ip-^ip-) and (■?/'+,'?/'+) fit into the same (2,2) supersymmetry 
multiplet as 0. 

2.1 The action 

The Euclidean action that we study is 

S = j (fz [-^d.d^cj) - 2ixp_dj^_ + 2i^+d,i)+ - FF 

(0)F + w\^)F - w"{(t))i)+i)- - W"(0)^_V'+j (2.1) 

In these expressions, z = xi + ix2 = (z)*, = (<9i — 182) /2 = (c^)*, d'^z = dzdz, F 
and (p are complex scalar fields, W'{(j)) = dW/dcp, etc., where W{(j)) is holomorphic 
in (f) and is the superpotential. and are independent Grassmann fields: the 
partition function is Z = J[(P(j) (PF (fip_ (fip+] exp(— S'), where d'^ip^ = dip- dip^, 
etc. 

Although the quadratic term of the auxiliary field F has the "wrong" sign, it can 
be formally integrated by analytic continuation. This allows for the elimination of F, 
leading to the action 

S = J d^z [-4:(j)d,ck(p - 2itp_ch^_ + 2i^+d,ij+ 

+ \W\^)\^ - W"{(P)ij+ij. - W"{^)ij.ij+j (2.2) 

The continuum partition function is thus Z = J[(f(p(f'ip_ c^^V'+] exp(— 5), now with 
action TTB . 

Of course Z is at this point only formally defined. Meaning can be assigned to it 
through perturbation theory about = 0, or some more general set of classical sad- 
dlepoints. In either case, one obtains an asymptotic series that provides approximate 
results valid in a limited regime. However it is worth noting that the supersymmetry 
of the model leads to nonrenormalization theorems that allow for a number of exact 
results to be obtained; see for example ^Hj and references therein. It is hoped that 
the lattice formulation that we describe below will give a more complete meaning to 
the functional integral, outside of perturbation theory and semiclassical expansion. 
In that case, aspects of the theory that are not protected by nonrenormalization 
theorems can be studied in a strongly-coupled or deep infrared regime. 

We speciahze to the superpotential 

W{<P) = + (2.3) 
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Both m and g have mass dimension 1. Perturbation theory is defined hj g/m 1. 
It is useful to make the field redefinition 



-^ + a (2.4) 
9 



In terms of a, the superpotential becomes 



W{a) = -\a + l^a^ A = ^ (2.5) 
2.2 Classical symmetries 

With periodic boundary conditions for all fields, the action ()2.1|) is invariant under 
the infinitesmal supersymmetry transformations (also true for (p ^ a) 

6ip. = 2ie'd,<p- e^F, SF = 2ie-d,ip+ + 2ie^chip_ 

5<j) = -e-^)^ - e+V'+, 5^'+ = 2ie^d^^ + e'F 
b^- = -2ie'd,<t)-e^F, 5F = 2ie-d,ip+ + 2ie^ch^_ (2.6) 

It can be seen that e~^,e~^ are associated with whereas e~,e~ are associated with 
dz. We can associate supercharge operators with the transformations ()2.fi|l through 

S = e-g_ + e+g+ - e-g_ - e+g+ (2.7) 

From this definition it is straightforward to work out the action of Q±, Q± on the 
fields. Two subalgebras emerge: 

{g_, g-} = -2tdz, {g+, g+} = 2^^ (2.8) 

All other anticommutators vanish. Note that from the relation tM = —'He between 
Minkowski and Euclidean time, 

z = X + He = X — tM, z = X — He = X + tM (2.9) 

Thus the first subalgebra in ()2.8j) closes on the generator of left-moving translation 
and the second subalgebra closes on the generator of right-moving translation. Re- 
spectively, there is a left-moving (2,0) algebra and a right-moving (0,2) algebra. The 
"2" denotes two generators, say g_ and (Q-- Taken together, the system has (2,2) 
sup er symmet ry. 

The superpotential ()2.5j) has a U{1)r symmetry when X = 0, g ^ (note that 
(j) = a when A = 0): 



a —>■ e 



ij± e-^°/V±, F e-^^^/^F (2.10) 
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and a,ip±,F transforming as conjugates. Note that a has R-charge 2/3. Thus W 
exp{2ia)W , which imphes that the superpotential (still with A = 0) has R-charge 2; 
this is just the well-known condition to have a U{1)r symmetry. That this is an R- 
symmetry follows from the fact that the component fields of the same supermultiplet 
have different R-charges, as can be seen from ()2.10|) . 

There is also an axial U{1)a symmetry that we will make use of: 

V^i^e^^-V'i, V-i^e^^^i (2.11) 

with all other fields neutral. Note that this contains fermion parity Z2{F) as a 
subgroup. 

We now show that the classical theory with X 0, g has an exact Z2{R) 
symmetry. The scalar potential is 

(2.12) 



\W'ia)\^ = 


A cr 
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V 

This potential has minima (recall A = m? /2g) 

a± = ±^2\/g = ±m/g (2.13) 

Clearly ()2.12|) has a Z2 symmetry a —>■ —a that relates the two classical vacua. The 
Z2 is a subgroup of the U{1)r x U{1)a that survives when A 7^ in ()2.5|) . as will be 
discussed below. Since V = \W'{a)\'^ > 0, these are absolute minima with V{(t±) = 0. 
The saddle point that separates these minima is at the origin, with 

I |4 

W"{0) = 0, V{0) = |A|2 = ^ (2.14) 

It can be seen that as \g\ is decreased with m held fixed, the minima separate and 
the energy density of the barrier between them increases in height. 

The effect of A 7^ is best understood by treating it as a background field with 
nonvanishing R-charge. For the superpotential to have the requisite R-charge of 2, A 
must have R-charge of 4/3: 

A ^ e^^^/^A (2.15) 

Now notice that A 7^ is left invariant iff 

37m 

a = neZ (2.16) 

Thus nonzero A leaves intact the Z^i^R) subgroup generated by a = 37r/2 in Eq. ()2.1U|) . 
At the level of component fields, this tranformation is (to = e™l'^\. 

o fi(c<j)cr = uP'a = —a, 4'± ~^ ^{uj)ip± = oj^^ip^ = —iip± 

F n{uj)F = F (2.17) 
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and a,il)±,F transforming as conjugates. Here, Q is the homomorphism that deter- 
mines how the Z^IK) acts on the fields; i.e., their 4-ahty. It is easily checked that 
this is a symmetry of the full action with superpotential (|2.5p . 

To obtain the promised Z2{R), we combine the transformation ()2.17|) with a Zi{A) 
subgroup of the U{1)a that appeared in in ()2.11|) . We describe it by the action of a 
generator 73: 

73 : iip+, ip- -i^J-, i>+ -ii>+, i>- i^- (2.18) 

The Z2{R) is generated by •y^Q^u): 

Z2{R) : a ^ -a, ^ -a, ^± ^ ±7/-^, ^± ^ ±i,± (2.19) 

Thus for A 7^ 0, (7 7^ 0, or equivalently m 7^ 0,(7 7^ 0, the surviving global symmetry 
is U{1)a X Z^{R) = U{1)a X Z^iR). 

At infinite volume, the classical Z2{R) symmetry must be spontaneously broken. 
This is because the scalar potential is positive semi-definite: V = |Vr'(cr)p. In this 
case, the stability of the scalar potential implies in 2d that there are no nontrivial 
finite action solutions: only the constant absolute minima a±. of ()2.13p have finite 
action.^ Thus instantons do not exist, tunneling does not occur, and the symmetry 
is spontaneously broken. 

2.3 Renormalization 

The only renormalization of the superpotential W is that due to wavefunction renor- 
malization jlj, with all component fields in the same supermultiplet rescaled iden- 
tically.^ Thus if the bare fields and renormalized fields are related by = \/Z(f)r, 
'ijj± = y/Z'ijj±^r, etc, then the renormalized superpotential is completely described by 
the identification W{m,g\(f)) = W{mr, gr\(pr) with 

TTlr = Zm, Qr = Z^^^g (2.20) 

Note that this implies = ^/ZX. Counterterms 5m and 6g for the couplings are 
not required. I.e., there is no additive renormalization of m or g; if the bare mass 
m = 0, then the theory remains massless under renormalization. This can also be 
seen as a consequence of the f/(l)_R symmetry ()2.10|) that occurs if m = 0. Since a 
continuous symmetry cannot be broken in 2d [TSJ, we know that a mass term cannot 
be radiatively generated, since that would break the U{1)r symmetry. Thus in the 
continuum theory the critical domain is the neighborhood of m = for any g. We 
will find that this is likewise true in the lattice theory. 

^See for example Section 41.4 of [T7| . 

^Of course, Wess and Zumino proved the nonrenormalization theorem in 4d; but the crucial 
cancellations, due to the 4 conserved supercharges, continue to hold in 2d. 
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3 The lattice action 



We now describe the Q-exact lattice action that is the subject of our Monte Carlo 
simulation studies. The general form of this action was first constructed in [311] using 
a Nicolai map [S] approach. In P, Monte Carlo simulations were performed. These 
demonstrated that Ward identities associated with the (2,2) supersymmetry of the 
continuum theory were satisfied to such a good approximation that no statistically 
significant violation could be observed, even with a large sample size. In |8j and 
[H], two different methods demonstrated boson-fermion spectrum degeneracy for the 
lightest states, a necessary condition for supersymmetry. Furthermore, the numerical 
values for the light spectrum agreed with those obtained on the Q, Q^-preserving 
spatial lattice in |3|. In |lTj the lattice action was derived using a lattice superfield 
approach. It was shown how to introduce a Wilson mass term in the superpotential 
to lift doublers, while keeping the one exact supersymmetry. The same Wilson mass 
term also appeared in the Nicolai map construction of O |3] . In ^T] it was shown that 
in the a —>■ limit the lattice perturbation series goes over to that of the continuum 
without the need for any counterterms, to all orders in the coupling g. 

With auxiliary fields eliminated, and implicit summation over repeated lattice site 
indices m, n, the action is:^ 

S = 4[Af0+^W^'(0)]„[A|0-^W'(0)]„-2#_,„A|^_,„ 

+2i?/'+,„AfV'+,„ - llj+^rnW"{(f))mn'ip-,n - ^ - ,niW" {(j)) mn^ + ,n (3.1) 

For convenience, we define the superpotential to include a sum over lattice sites: 

^ = Ef-i<^-^'<^- + E|c) (3.2) 

The purpose of the Wilson mass term </)A^</), which includes nearest neighbor inter- 
actions, is to lift doublers in the spectrum in a supersymmetric way. From ()3.2|) one 
obtains 

dW „ d'^W 

dcpra 0<pmO<pn 

In expressions fj3.1|) and ()3.2p . the following finite difference operators are used: 

= 1 f A+ + A") , A^ = V A+A- 

fi=l,2 

Af = ^ (Af - ^Af ) , A| = i (Af + ^Af ) (3.4) 



*In all of our expressions we work in the natural lattice units a = 1, unless otherwise noted. 
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where and A~ are the usual forward and backward difference operators respec- 
tively. From now on we drop the superscript S, leaving it implied, except on the 
operator 

A^^ = i:^'^' = 4AfA| (3.5) 

which we want to distinguish from A^. We will often suppress site indices as well. 

We pause to note that ()3.1|) may be obtained in terms of the one exact lattice 
supersymmetry Q acting on an expression of component fields. This Q is the linear 
combination Q = Q_ + Q+ of continuum supersymmetries ()2.(j|l - ()2.7j) . with dz — > A^ 
and dz —>■ A^. Note that 

Q' = l{Q- + Q+,Q- + Q+} = (3.6) 

so that Q is indeed nilpotent. For reference, the action of Q on the component fields 
is 

Q(l) = tlj-, = F + 2iAz(j), Q^p-=0, QF = -2iAzij-, 

Q4) = = 0, Q^. = F- 2zA,0, QF = 2iAzi>+ (3.7) 

from which one can explicitly verify the property (j3.6j) . For a general holomorphic 
function W^(0) we find that 

S = g(-F^_ -2zV^+A,0 + l^'(0)7/^++W'(0)7/;_) 

= -2#_ AjV'- + 2iilj+Azip+ - 4:(j)AzAz(f) -FF - ij_w"i}+ 

-ip+W"^. + W'{(j)) {F + 2iAz(j)) + W(0) {F - 2iAz^) (3.8) 

Upon elimination of the auxiliary field F through its equation of motion, (j3.8|) be- 
comes just the action p.ip . 

Because ()3.8|) is Q-exact, i.e., S = Q{- ■ •), we know that QS = 0. That is, the 
exact supersymmetric invariance of the lattice action w.r.t. Q just follows from the 
nilpotency of Q and the Q-exactness of the lattice action. This obviously provides 
a mechanism by which a host of supersymmetric lattice systems can be constructed. 
In pi] it was shown that (j3.1|) is far from the most general action with the symmetries 
of 1)3.11) . Hence the importance of studying the renormalization of the lattice theory. 
In particular, we would like to examine whether or not relations like ()2.20|) hold. In 
our simulation results below, we will find evidence that ()2.20p is approximately true, 
in a regime where a <C 1. 

It is a simple matter to specialize to the case of the cubic superpotential. Then 

^('^) = E {-l4>n.A'4>^ + + (3.9) 
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The action can be written as = 5*0 + Sint-, where 5*0 is quadratic in fields: 



Sn 



-A^2 + (m-^A2)(m-^A2) 



+2#+A,V^+ - ^+(m - 2^')^- - ^-("i - (3.10) 

At r = the bosons and fermions both have doublers. The Wilson mass terms at r 7^ 
lifts these modes. Note that the inverse free propagator for the boson is the square of 
that for the fermions. This is a consequence of the exact lattice supersymmetry, and 
is a property that is crucial to the perturbative results that were mentioned above. 
The interaction terms are contained in 

- - - 1 - 1 -1 

Sint = -g4>i'+^P- - g(pi^-i'+ + -gm4)^(p + -9rh4>^(p + \g?(t>^4>'^ 

-^r#2A2^ _ lfg4>^^^ + ig<f^^<P - «^0'A,0 (3.11) 

Note that irrelevant operators appear in the second line. These are a consequence of 
the exact lattice supersymmetry. The r-dependent terms originate from the super- 
symmetrization of the Wilson mass term. The other two arise from crossterms in the 
first line of ()3.H) . such as Vr'(0)A20. The crossterms only contribute to interactions, 
due to the identity 0A^0 = on a periodic lattice. 

The irrelevant operators coming from crossterms violate the rotation symmetry 
of the lattice, as well as refiection positivity. As is well-known, Z4 rotation symmetry 
is sufficient to imply Euclidean rotation invariance, 5*0(2), of the continuum limit. 
Thus the violation of Z4 rotation symmetry is an important matter, since it could 
conceivably destroy Euclidean invariance of the continuum limit. It can be shown JT] 
that the breaking in this system is an 0{a?gphys) effect, where g^hys is the coupling 
in physical units. It is harmless to the continuum limit in perturbation theory, as the 
lattice perturbation series turns out to be finite. However, it remains to be shown 
that this does not affect the continuum limit at a nonperturbative level. We will defer 
detailed discussion of this issue to a later publication. Suffice it to say, we have found 
in simulations that Green function data at more than a few lattice spacings exhibits 
approximate Z4 symmetry, and that this approximation improves as one approaches 
the continuum limit. As to the lack of refiection positivity, this too is an 0{a^gphys) 
effect. Some problems were detected in propagators at gphys = in JS]; however, 
these difficulties apparently went away for sufficiently small lattice spacing a. In the 
studies that we perform here, we work in a regime where g = gphysO- <^ 1, where the 
violation of refiection positivity by irrelevant operators is believed to be harmless. 

Next, consider the field redefinition (j2.4|) in the lattice theory. One finds 

W{a) = (-^a^AV™ - Xam + (3.12) 

m 

and the lattice action is just ()3.1|1 with (p a. If r = 0, then W'{a) is even 
and W"{a) is odd under the Z2{R) transformation ()2.19j) . Taking into account the 
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transformation of the fermions, the action is by inspection Z2{R) invariant. However, 
the Wilson mass term violates the Z2{R) symmetry (|2.19|) since for r 7^ 0, W'{(j) and 
W"{a) no longer have definite parity with respect to Z2{R). On the other hand, the 
symmetry breaking arises entirely from irrelevant higher-derivative terms. Given that 
the lattice perturbation series is finite, the Z2{R) is restored in the continuum limit 
in perturbation theory. What is left to check is that this occurs nonperturbatively. 
We will find rather strong evidence that this is true. 

Finally, we note that in all of our simulations we set r = 1 in the Wilson mass 
term. 

4 Simulation results 

4.1 Location of the critical domain 

First we introduce observables that allow us to locate the critical domain.^ For 
a generic choice of the bare parameters m and g, the effective potential has two 
minima, in correspondence with the classical minima o'± of ()2.13|) . We will refer to 
these two minima of the effective potential as the a±. vacua, although their location 
will be somewhat shifted from the classical values. The definitions that we choose 
are intended to address the following issue related to the existence of two vacua. At 
finite volume tunneling between the two vacua a± will occur. This is important, 
because the tunneling events are very large fluctuations, of order m/g. We do not 
want to include such fluctuations in the measurement of correlations, since they do 
not correspond to fluctuations about one of the Z2{R) symmetry breaking vacua in 
the thermodynamic limit. By a± domains we will mean connected regions where a is 
closer to one minimum of the effective potential than the other. Several a± domains 
may establish themselves on the N x N lattice, particularly for large A^. In that case 
what we really want to measure is the fluctuation away from the average value within 
a domain. 

Under the assumption of approximate Z2{R) symmetry of the effective potential 
at finite volume, we expect that the following occurs. Take m, g real so that the 
classical minima ()2.13p are real. Denote the real and imaginary parts of a-m as aR^m 
and aj^m respectively. Suppose we have a ensemble of configurations F. We partition 
this ensemble as F = F_|_ + F_ + Fo where: aR^m > on the subensemble F_|., a^^rn < 
on the subensemble F_, and cTR^rn = on the subensemble Fq. Note that Fq is a set 
of measure zero in any expectation value over F, since (Tpt^m is a continuous variable. 
Let {o'r)_^_ be the average of an^m over the subensemble F+ and {<Jr)_ be the average 
for the subensemble F__. The order parameter (IctrI) will denote the average of |cr/j_m| 

'"^Of course in our analysis we work at finite volume, so all so-called critical behavior is actually 
pseudocritical. We will not belabor this point in our terminology by appending "pseudo-" to every 
instance of "critical." 
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N {aR)^ M_ {WrD 

4 3.505(52) -3.518(52) 3.512(14) 

8 3.337(49) -3.353(48) 3.3456(87) 

16 3.342(60) -3.349(41) 3.3471(45) 



Table 1: Comparison of subensemble order parameters to that of the full ensemble, 
for (m,c/) = (0.10,0.03). 

over the full ensemble F. Then we expect 

{aR)^ = -{an)^ = {\cTR\)^v (4.1) 

to a very good approximation. The veracity of this conjecture is borne out in our 
simulations. As an example, in Tabled we show results for the quantities of ()4.H1 for 
{m,g) = (0.10,0.03) and three different sizes of N x N lattice. Similar behavior was 
observed at other values of {m,g,N). This is already good evidence that Z2{R) is 
a symmetry of the effective potential, to within the 1-2% statistical uncertainties of 
the measurements. 

This motivates the following fluctuation analysis. In either a± domain (connected 
regions where cr > or cr < 0) we define SaR^m by \<7R^m\ = v + SaR^m- Thus the 
fluctuation 6aR^rn is related to the value of cr/j,m through 

(^R,m = V + 6aR^rn, R,m > 

0-R,m = -V - SaR^rn, (yR,m < (4.2) 

In either case, the fluctuation is positive if it moves in the direction that is away from 
the origin aR^m = 0, relative to ±v; i.e., in the direction that would increase |o"ij,m|. 
This is shown in Fig. [H In the thermodynamic limit, l^aij^ml <^ v and SaR^m is just 
a fluctuation about one of the two vacua, which is what we want to study. 

Associated with the fluctuation 6aR^rn one has the following (momentum space) 
Green function and susceptibility: 

GkiSaR) = ^{\6aR,k\'), xi^aR) = Go((5a«) (4.3) 

Here, V = N1N2 on an Ni x N2 lattice (typically we take Ni = N2 = N) and our 
convention for the Fourier transform is 

T~ e o .fkimi k2m2\ , , 

ocrR,k = OaR^rn exp 2711 I -j^ + -j^ 1 (4.4) 

A rough estimate of the correlation length is the one used in JHl- It involves the ratio 
of the Green function Gk{SaR) at k = (0, 0) to k = (1, 0): 



G{o,o)iScrR) 
4sin2(7r/Ari) \ G(i,o)(5a«) 



1/^ 



- 1 



(4.5) 
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Figure 1: Positive fluctuations San^m away from vacua ±v. 

Near the critical point we expect this quantity to exhibit the correct scahng behavior, 
at leading order. Since we are not interested in high precision spectral information, 
()4.5|) will suffice for our purpose (determining the location of the critical domain). 

We also want to consider fluctuations from the imaginary part of the fleld a. 
Numerically we have found that {(Ji^m) = 0, for any real m,g. We therefore deflne 
Sai-m. = (Ti^rn and 

5am = SaR^rn + «5(T/,m (4.6) 

This leads to the Green function Gk{Sa), susceptibility xi^^^) and correlation length 
C{Sa), obtained through the replacement 6aR 6a in ()4.3|) . ()4.5|) above. 

In Tables 12111] we display the measured value for the observables at g = 0.03, as 
a function of m, on three different N x N lattices, = 4, 8, 16. A value <^ 1 has 
been chosen because, as the reader will recall, g is the coupling in lattice units. Thus 
in terms of the coupling gphys in physical units, g = gphysO^- For gphys = 1, the choice 
g = 0.03 corresponds to a = 0.03. Similarly, m = rriphysCi, from which one obtains a 
physical interpretation of the mass. 

Several things are to be noted from the data: 

• The measured value of {\aR\) is close to the classical value \a±\ = m/g in the 
case g <^ m, where perturbation theory applies. Thus the sample is dominated 
by conflgurations for which am sits near one of the two classical minima, as one 
would expect. The classical estimate breaks down as m ^ 0. Note, however, 
that if we follow the data at m = g for = 4,8,16, the value of (|cr/j|) is 
decreasing, and for A^ = 16 is only 41% larger than the classical prediction. This 
implies that for sufficiently large A^, perturbation theory should be trustworthy 
in the regime g/m ^ 1, as one might guess. 

• The 6aR observables of susceptibility and correlation length indicate that the 
critical domain is in the neigborhood of m = 0. There is a plateau that is fairly 
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m 


{WrD 










0.00 


2.360(30) 


7.93(48) 


28.8(3.4) 


12.51(21) 


150.7(4.3) 


0.03 


2.462(14) 


7.86(22) 


29.5(1.6) 


12.060(93) 


139.8(1.9) 


0.05 


2.624(14) 


7.91(23) 


29.8(1.7) 


11.464(96) 


128.1(2.0) 


0.10 


3.512(14) 


7.91(28) 


31.2(2.2) 


9.23(13) 


84.8(2.3) 


0.15 


5.012(22) 


6.68(68) 


23.7(4.8) 


6.96(35) 


48.5(4.8) 


0.20 


6.6473(89) 


4.85(52) 


12.8(2.7) 


4.89(26) 


25.3(2.7) 



Table 2: Observables on the N — A lattice at gi = 0.03. 

(WrI) ^{^(^r) xiS(TB) ^(Sa) xiSa) 



0.00 
0.03 
0.05 
0.07 
0.10 
0.15 



1.593(29) 
1.8162(94) 
2.0852(99) 
2.5037(98) 
3.3456(87) 
5.0026(60) 



11.6(1.0) 
11.08(32) 
11.53(36) 

11.17(41) 

9.72(53) 

7.21(88) 



54.3(9.1) 
57.0(3.2) 
62.3(3.8) 
61.4(4.4) 
48.2(5.1) 
22.7(5.3) 



18.30(41) 
16.68(14) 
15.13(15) 
13.17(18) 
10.06(25) 
7.17(44) 



306(12) 
262.6(3.8) 
221.0(4.2) 
170.6(4.6) 
103.3(5.2) 
45.2(5.4) 



Table 3: Observables on the N = 8 lattice at g = 0.03. 



m 








((()a) 


a(^^^^) 


0.00 


1.2261(62) 


15.67(49) 


98.3(5.9) 


25.37(19) 


573.7(7.4) 


0.01 


1.2309(89) 


15.93(70) 


100.9(8.5) 


25.29(27) 


564(11) 


0.02 


1.2802(89) 


15.86(71) 


102.2(8.8) 


24.21(27) 


526(11) 


0.03 


1.4076(66) 


16.16(53) 


110.9(7.0) 


22.73(20) 


472.5(8.0) 


0.05 


1.7818(96) 


15.88(86) 


116(12) 


18.60(37) 


330(13) 


0.07 


2.3579(86) 


13.9(1.1) 


95(14) 


14.56(52) 


208(15) 


0.10 


3.3471(45) 


9.9(1.1) 


51(11) 


9.89(56) 


101(11) 



Table 4: Observables on the N — lattice at 51 = 0.03. 
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flat (within statistical uncertainties) extending out from m = to some value 
of m that is not too large. The range of this plateau appears to be shrinking as 
is increased. From this we conclude that the critical domain shrinks down 
to an infinitesmal neighborhood of m = in the thermodynamic limit. Note 
that this is consistent with the nonrenormalization theorem of the continuum 
theory, expressed by ()2.20j) . That is, if the bare mass m in fl2.20|) is set to zero, 
the renormalized mass is also zero: there is no additive renormalization. 
As mentioned above, this can also be viewed as a consequence of the f/(l)R 
symmetry, which protects against a mass term in the superpotential. The fact 
that on the lattice the critical domain is in a shrinking neighborhood of m = 
leads us to expect that the U{1)r symmetry holds to a good approximation for 
the IR modes in the continuum limit. Further evidence for this will be seen 
below in Section lOl 

• There is a hint of a dip in xi^'^fi) at m = 0, but the statistical errors are too large 
to be certain. Simulations at m = are very costly, since the leapfrog evolutions 
in the hybrid Monte Carlo algorithm become quite unstable, as discussed in 
Appendix |X] In any case, the definition ()4.2|) of da^ does not make much sense 
right at m = 0, since (at least classically) there is no barrier and the vacuum is 
unique. At m = it would be more appropriate to study fluctuations of (jR^m 
rather than |o"ij^m|- However, since our purpose here is only to find the critical 
domain, it is enough to see that the susceptibility x{^'^r)j which does make 
sense at m 7^ 0, is rising as we approach m ^ 0. 

• The Sa observables do not show a plateau, but are maximized at m = 0. The 
fluctuations Saj contribute significantly, and in fact dominate, as follows from 
x((5(t) = x{^(^r) + x(^cr-f)- The degree to which this is true increases as we 
approach m = 0. We will find further evidence below that m = is the critical 
point. As a result, the 6a observables seem to be better indicators of criticality. 

In conclusion, the critical domain is the neighborhood of m = 0. The continuum 
renormalization ()2.20|) appears to hold to a good approximation. At m g, the 
theory is noncritical. 

4.2 R-symmetry of the effective potential 

Next we address the two principal questions of this work: (1) To what extent is 
Z2{R) a symmetry of the effective potential? (2) Is ?7(1)/? symmetry approximately 
recovered in the critical domain? 

We introduce a complex constant external field h 




(4.7) 



m 
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to the scalar potential V. This allows us to explore the effective potential and to study 
the extent to which it is symmetric w.r.t. cr — > — cr, or the phase rotation a e^^a. 
As usual, we obtain the effective potential from the Legendre transformation of the 
generating function w{h) = In Z[h), where Z{h) is the partition function that is 
obtained when ()4.7p is added to the lattice action. Z2{R) symmetry of the effective 
potential is equivalent to w{—h) = w{h). Similarly, U{\)r symmetry of the effective 
potential is equivalent to w{e^^h) = w{h). Note also that (cr)^ = dw{h)/dh, where 
(o")^ is the expectation value of a in the background h. It follows that in the case of 
Z2{R) symmetry we have the prediction 

{^)-h = -{^)h (4.8) 
In the case of symmetry we have the much stronger prediction 

{a)^,e, = e^'{a), (4.9) 

Equivalently, since we take m > 0, g > and will find below that {a)f^ > if h is real 
and positive, 

arg(cr)^ = axgh, |(cr)^| = const., fixed \h\ (4-10) 

is the prediction of f/(l)/j symmetry. The second condition just states that Ko")/^^ = 
(cT/j)^ + {'^i)'h independent of arg/i. 

4.2.1 Z2{R) tests 

Here we consider the case where h is real. It can be seen from Table El that, up to 
statistical errors, the values of {an) are supportive of the Z2{R) symmetry prediction 
(j4.8|) . It is a bit surprising that this is already true on such a small lattice [N = 4). Of 
course g/m = 0.15 is well within the perturbative regime, so upon further reflection 
it is reasonable that the irrelevant Wilson mass interaction terms that break the 
Z2{R) symmetry would only have a mild influence on the effective potential. Still, 
it is remarkable that at finite lattice spacing the irrelevant operators have negligible 
effect. 

In Fig. |21we show the behavior of {cr),^ on lattices with {m,g) = (0.20,0.03) and 
= 4, 8, 16. It can be seen that the transition sharpens significantly as we go from 
the = 4 lattice to the A^ = 16 lattice. Due to nonergodicity in the simulations 
when the barrier between vacua is large, it was not possible to get stable averages 
in the transition region, very close to h = 0, for the A^ = 16 lattice. So, we have 
performed a hysteresis study at {m,g,N) = (0.20,0.03, 16). 

A simulation was done where a thermalized configuration with h = —0.1 was 
first produced. Then (aR) was computed from successive configurations obtain in 20 
updates. Next, h was increased by Ah = 0.01. Then 20 updates were performed 
and (ctr) was calculated from the 20 configurations that were generated. This was 
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h 



{aR){h<0) {aB){h>Q) 



0.10 
0.01 
0.005 
0.003 
0.001 
0.0005 
0.0001 
0.0 



-8.427(19) 8.405(17) 

-6.663(61) 6.7231(45) 

-5.395(56) 5.428(53) 

-3.776(79) 3.897(79) 

-1.408(95) 1.520(98) 

-0.667(99) 0.694(25) 

-0.235(96) 0.145(96) 

-0.038(97) — 



Table 5: {crR)i^ as a function of the source h for (m, g, N) = (0.20, 0.03, 4). Note that 
{(Ji)^ = for /i reaL 

continued up to /i = 0.1. Then the cycle was reversed, decreasing h by the same incre- 
ment. The results are presented in Fig. El It can be seen that there is a pronounced 
hysteresis. We increased the number of successive configurations at each step of the 
cycle from 20 to 200 and found that the hysteresis loop did not tighten at all. This 
is indicative that the transition at ^ is first order. This is as we should expect 
from Table m the point {m,g) = (0.20,0.03) is well outside the critical domain. It is 
only in the critical domain that the two vacua begin to degenerate, giving rise to the 
critical fluctuations which occur in a second order transition. 

Thus, an interesting question is the extent to which the first order transition 
softens in the vicinity of the critical domain m ^ 0. To study this we have produced 
hysteresis curves like Fig. Elfor a sequence of decreasing m, holding g = 0.03 fixed. It 
can be seen that the result is entirely consistent with the conjecture that the critical 
domain occurs in the neighborhood of m = 0. The hysteresis curves close up as m is 
reduced, indicative of the onset of a second order transition. 

Note that the h = data. Ah > versus Ah < 0, are closer for m = than 
for m = 0.05. This supports the conclusion that m = is the critical point, or the 
center of the critical domain. From this we draw a conclusion about the plateau and 
possible dip in 6aji observables that was noted above in relation to Table El versus 
the m = maximum in 5a observables, also seen in Table 01 The fact that less 
hysteresis is observed at m = than at m = 0.05 indicates that the former point is 
"more critical." Since the 6a observables are noticeably peaked at this point, whereas 
the San observables are not, we conclude once again that the former observables are 
better indicators of criticality. Yet more support for these conclusions will be obtained 
below. 
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Figure 2: (ct/j) versus h for {m,g) 



(.20, .03). 



4.2.2 U{1)r tests 

Our tests of the t/(l)_R conjecture ()4.10|) were conducted as follows. Very large en- 
sembles (10"^ independent configurations) were generated at {g,h,N) = (0.03,0,16), 
for various values of m. Expectation values {o')f^ at nonzero h were then obtained by 
reweighting IQ^ with the potential ()4.7|1 : 



The accuracy of the reweighting procedure was verified by comparison to simulations 
run at a nonzero values of h. Note that the reweighting (j4.1ip was not used in the 
Z2{R) tests of Section 14.2.11 above. There, simulations included the potential (jUTj). 
Comparing to expectation values obtained there gives us confidence in the reweighting 
procedure 1)4.111) . 

We keep h small for three reasons. First, to maintain the reliability of the reweight- 
ing procedure, since a large value of h will tend to lead to overlap problems; i.e., the 
ensemble of configurations simulated at h = may be distributed very differently 
from one with a large value of h. Second, for the = 16 lattice where we will 
test for U{l)ii, the volume A^^ amplifies the effect of h on the homogeneous mode, 
ak=o = Ylim^rn) whosc bchavior is characterized by the effective potential. There is 
already a large sensitivity in (cr)^ = N~'^{ak=o) at very small h. Third, what we are 
after are the symmetries of the effective potential in the h ^ limit. 

In Fig.Elwe display a.Tg{a)i^ versus arg/i at {g, |^|,iV) = (0.03,0.001, 16) for three 
different mass values, m = 0,0.03,0.10. For m = 0, the data passes through the 



(4.11) 
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Figure 3: Hysteresis cycles for (aji) versus h at {m,g,N) = (0.20,0.03, 16). The top 
figure had 20 successive configurations at each value of h along the cycle, whereas the 
bottom figure had 200. 
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Figure 4: Hysteresis cycles for (aji) versus h at {g, N) = (0.03, 16), for various values 
of m. We averaged over 20 successive configurations at each value of h along the 
cycle. 
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(diagonal) straight line arg (a) = arg/i, showing that f/(l)_R is a very good symmetry 
of the effective potential. For m = 0.03, the data deviates slightly from the straight 
line, indicating that the symmetry is only slightly violated. Finally, at m = 0.10, 
the symmetry is completely broken. The fact that arg (a)^ ~ ±7r in this case can be 
understood as follows. For larger values of m and the very small h that we choose, the 
potential V = \W'{a)\'^ dominates over the source potential AV{h) of ()4.7p . In that 
case, (cr)^ ^ ±v. Recall that v was defined by (j4.1|) . without the source potential. 
Also recall that v is real for m, g real, as we choose here. The role of h then is just 
as a perturbation to pick the sign of ±v. It follows that arg (cr)^ ~ ±n. 

The correlation between Fig. El and the results of Table 0] is very good: on the 
plateau where Sa^ observables take their critical values, the U{1)r symmetry is a good 
approximation. At the point where the 6a observables are maximized, no f/(l)i? 
symmetry violation can be seen in Fig. El above the statistical uncertainty. This 
indicates that m = is indeed the critical point, and shows that 6a observables are 
well-suited to studying behavior at that point. We again find that any possible dip in 
the 6a R observables at m = has no significance to the location of the critical point. 

The second part of the conjecture (j4.10|) was studied through the quantity 

Rmj) = ^-^^^, M = ^EK-).,i (4.12) 

where hj = \h\ exp(27rij/n) corresponds to the values of h that were used in the data 
set. Thus, R measures the relative shift of away from the mean, where the mean 
is taken w.r.t. arg/i. In Fig.lHlit can be seen that significant deviations away from the 
mean are observed at m > 0, with the largest occuring for the value m = 0.10. (The 
R test of Fig. ini turns out to be more sensitive to the U{1)r violations at m = 0.03 
than the arg (a)^ vs. arg h test of Fig. Elwas.) On the other hand, it can be seen that 
the U{1)r symmetry conjecture (j4.10|) holds up quite well at m = 0. Once again, 
m = is found to be "more critical" than m ^ 0.03, and the 6a observables seem to 
be better indicators of criticality than the 6aR observables. 

In conclusion, the U{1)r symmetry of the effective potential is quite robust at 
m = 0. To within statistical errors it is not violated. This provides strong evidence 
that additive renormalization of mass does not occur in the lattice theory in the 
continuum limit. The violation of the ?7(1)h symmetry in the lattice action at m = 
due to the Wilson mass term in the superpotential does not appear to impact the 
renormalization properties of long distance modes. 

5 Interpretation and future directions 

The simulation results that were presented above are quite encouraging for the Q- 
exact action of the (2,2) 2d Wess-Zumino model. We have found that the R-symmetry 
of the continuum theory is a property of the effective potential in the lattice theory. 
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Figure 5: A test of symmetry, by comparison to the prediction arg (cr)^ = arg h, 

indicated by the (diagonal) straight hne (mostly hidden by the m = data). It can 
be seen that the m = data is in very good agreement with the conjecture of U{l)ii 
symmetry. At m = 0.03, where susceptibilities and correlation lengths are close 
to their critical values, the f/(l)_R symmetry is approximate. At m = 0.10, where 
from Table m we see that we are certainly outside of the critical domain, the f/(l)R 
symmetry is badly broken. This data is for {g, \h\, N) = (0.03, 0.001, 16). 
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Figure 6: A test of t/(l)_R symmetry, by comparison to the prediction that \{cr),^\ 
should be independent of arg/i. Relative deviation from the average w.r.t. aigh is 
measured by R (c.f. eq. ()4.12p ). The prediction of U{1)r symmetry is indicated by 
the straight line R = 0. Again, the m = data is in very good agreement with 
the conjecture of symmetry, m = 0.03 data indicates that the symmetry is 

approximate (a maximum of 25% relative violations), whereas m = 0.10 data shows 
the symmetry is badly broken. This data is for {g, \h\, N) = (0.03, 0.001, 16). 
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The explicit breaking of the R-symmetry due to the Wilson mass term in the super- 
potential is harmless in the continuum limit; the continuum R-symmetry is recovered 
without the need for counterterms. 

A very important consequence of this is that the nonrenormalization theorems of 
the continuum theory appear to hold at a nonperturbative level for the long distance 
modes of the lattice theory. Thus, the good behavior that was already seen to all 
orders in perturbation theory ^Tj seems to persist in the strongly coupled regime 
g/m ~ 1, where renormalizations are significant. Together with the supersymmetry 
Ward identity results of j9j and the spectral results of |8,^, the results presented here 
provide strong evidence that the nonperturbative physics of the continuum theory 
can be reliably studied with the Q-exact action. This, in spite of the fact that it 
only preserves one out of four supercharges and explicitly breaks the R-symmetry, 
Euclidean invariance, and reflection positivity that are so important to the continuum 
theory. 

Undoubtedly these positive results are related to (i) the fact that the symmetry 
breaking is due to irrelevant operators and (ii) that IPI diagrams of UV degree 
D > do not occur in the lattice perturbation series. The cancellations of D = 
contributions of subdiagrams in lattice perturbation theory is intimately related to 
the exact lattice supersymmetry It would be very interesting to know whether or 
not other lattice actions with an exact supersymmetry, such as the super- Yang-Mills 
examples that have been recently proposed ^31201; have a finite lattice perturbation 
series, in the sense that they have no D > IPI diagrams. However, a careful power- 
counting analysis, comparable to that done by Reisz for 4d Yang-Mills j2I], has yet to 
be performed.^ If in these lattice super- Yang-Mills theories the perturbation series is 
finite, the results presented here give some hope that those theories have the correct 
continuum limit at the nonperturbative level as well. 

In research that is in progress, we are currently subjecting the Q-exact lattice 
system studied here to another nonperturbative test of its continuum limit. As men- 
tioned at the outset, the continuum theory in the critical domain is believed to afford 
a Landau- Ginzburg description of the minimal discrete series of A/" = 2 superconfor- 
mal field theories. As a result, the critical exponents of all relevant operators are 
known exactly. If the lattice theory has the correct continuum limit, it should be 
able to reproduce these exponents. We are presently studying this issue through the 
examination of hyperscaling (dependence on correlation length) and finite-size scaling 
(dependence on system size) behavior in the critical regime. We hope to report the 
results of that study in the near future. 



^As is well-known, the UV behavior of the lattice perturbation series is often worse than the 
continuum series, due to modified vertices and propagators. 
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Appendix 

A Simulation details 

In our simulations, we want to study the behavior of the theory in the critical domain. 
However, in that case Monte Carlo simulations tend to exhibit critical slowing-down 
due to long wavelength modes that do not efficiently decorrelate. The characteristic 
autocorrelation time t^c in simulations grows like 

TAc--\m-m,\~'^\g-g,\~'', z^,Zg = 0{l) (A.l) 

as we approach a critical point {m^,g^). This is a consequence of the local trajectories 
through configuration space that a typical simulation method uses. In some systems, 
global moves [22^ based on percolation algorithms |^ can be engineered to completely 
overcome this problem. For example, in scalar field theories with a Z2 symmetry 
one can exploit cluster algorithms that have been developed for the Ising model 
|22[ 1^ . Such a global method is not known for the present system.*^ Instead, we 
have used the method of Fourier acceleration [2^1, as applied fSEl EI! to the hybrid 
Monte Carlo (HMC) algorithm |2H]- With Fourier acceleration, different Fourier 
modes are evolved with disproportionate time steps during the leapfrog trajectories 
of the HMC algorithm. This allows the longest wavelength modes to be displaced 
farther through configuration space, speeding up their decorrelation. Simultaneously, 
shorter wavelength modes are displaced less far. This keeps acceptance rates in the 
Metropolis step of HMC high. 

We have found that Fourier acceleration significantly reduces autocorrelation times 
for the present system, as previously noted by Catterall et al. 0. At points where 
we have compared the performance of conventional leapfrog to Fourier accelerated 
leapfrog, we find a reduction of autocorrelation times by an order of magnitude or 
more. Furthermore, with Fourier acceleration applied to our simulations, we did not 

''It is possible that one might be able to exploit the approximate Z2 {R) symmetry of the effective 
potential that is present in this lattice model. A percolation method based on this Z2{R) might 
work well as a suggestor algorithm for Metropolis steps. I.e., a percolation method could be used to 
suggest new configurations, which would be accepted or rejected with a Metropolis criterion. 
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encounter problems with critical slowing-down. (However, see the discussion below 
regarding leapfrog integrator instabilities in the critical domain.) 

In Eq. (D.16) of [TTlj it was shown how to write the fermion matrix in a basis where 
it is real. Once this is done, we can introduce real pseudofermions yi^m, « = 1,2 to 
reproduce the fermion determinant: 



The pseudofermion representation permits HMC simulation of the system. We intro- 
duce momenta p and vr, conjugate to and y resp., and an auxiliary Hamiltonian: 



With Fourier acceleration, the HMC algorithm invokes molecular dynamics trajecto- 
ries that evolve the system according to a modified version of Hamilton's equations. 
The net effect is that the leapfrog steps are done in Fourier space, with each mode 
evolved with its own molecular dynamics time step. For further details, see the dis- 
cussion in |77j . 

After the random update of the conjugate momenta that occurs prior to each 
molecular dynamics trajectory, we evolve the Fourier modes, labeled by k, for T time 
steps of spacing dtk- Thus there is a molecular dynamics evolution for simulation time 
Tk = T ■ dtk between each randomization of the momenta. The Fourier acceleration 
occurs by choosing dtk oc l/tUfc, where Uk is determined by the frequency of modes 
in the — > limit, with some effective accelerator mass rriacc, as described in j27j . 
Although the random updates of the momenta introduce noise, in a generalization 
of the discrete Langevin equation (see for instance the very nice discussion in |,29j), 
it is still possible that some modes may not decorrelate (over reasonable periods of 
simulation time). This is due to the coincidences pointed out by Mackenzie ISO]- To 
overcome this potential problem, dt was randomized. We effect this in the following 
way. For a given number of steps T, we choose dtmax such that T ■ dt^ax = 7r/2. 
We then randomize dt uniformly in {0, dtmax)- The time step for each mode is then 
determined by dtk = dt/ujk- This time step is used for the entire molecular dynamics 
trajectory of T steps. This algorithm is motivated by the observation j^H] that in the 
(7 = theory it can be shown that the autocorrelation of a mode k is proportional 
to cosT{TdtkUJk), where cosy is an approximation to cosine that becomes exact as 
T — i> oo. Thus by choosing TdtkUJk = 7r/2 we get almost complete decorrelation in 
the g/m <^ 1 regime. In our simulations we find rapid decorrelation outside of this 
regime as well. 

We have measured autocorrelation times throughout all of our simulations. We 
find that 



Sp = \y^{M^M)-^y 





(A.3) 



TAC = Oil) X T 



(A.4) 
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That is, a few hybrid molecular dynamics trajectories are sufficient to decorrelate 
the configurations. (Hybrid refers to the randomization of momenta prior to each 
trajectory.) Non-ergodicitites were excluded by verifying that measured quantities 
were independent of initial conditions, and that subensembles gave identical results, 
within statistical uncertainties. 

Unfortunately, as we approach the critical point m = 0, the number of steps T 
in the leapfrog evolution must typically be chosen very large. Recall that Tdtmax 
is held fixed. If dtmax is too large, we find that the leapfrog integrator can become 
unstable jHHES- Thus in practice we choose dtmax as large as leapfrog stability will 
allow, and adjust T such that Tdtmax = 7r/2. In the critical domain, acceptance 
rates are then found to be nearly 100 percent. (Therefore the only purpose of the 
Metropolis step of HMC is to keep the algorithm "perfect.") In practice we find that 
at or near m = 0, it is necessary to have T = (9(10^). Thus if simulation time is 
measured in units of dtmax ^ we do not escape critical slowing down; this is because 
dtmax niust be drastically reduced as the degree of criticality (measured, say, by the 
correlation length in units of lattice spacing) is increased. A significant speed-up 
in the simulations could presumably be achieved through the multi-pseudofermion 
methods discussed in [221 ■ 
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